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Abstract 

We study the quantized Friedmann-Lemaitre-Robertson- Walker (FLRW) model minimally cou- 

n 

pled to a free massless scalar field. In a previous paper, solutions of this model were constructed 
as gaussian superpositions of negative and positive modes solutions of the Wheeler-DeWitt equa- 
tion, and quantum bohmian trajectories were obtained in the framework of the Bohm-de Broglie 
(BdB) interpretation of quantum cosmology. In the present work, we analyze the quantum bohmian 
trajectories of a different class of gaussian packets. We are able to show that this new class generates 
bohmian trajectories which begin classical (with decelerated expansion), undergo an accelerated 
expansion in the middle of its evolution due to the presence of quantum cosmological effects in this 
period, and return to its classical decelerated expansion in the far future. We also show that the 
relation between luminosity distance and redshift in the quantum cosmological model can be made 
close to the corresponding relation coming from the classical model suplemented by a cosmological 
constant, for z < 1. These results suggest the posibility of interpreting the present observations of 
high redshift supernovae as the manifestation of a quantum cosmological effect. 
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I. INTRODUCTION 

Recent measurements of high redshift supernovae 

aa 

indicate that the Universe is 

presently in accelerated expansion and not in a decelerated one, as was firmly believed by 
cosmologists since Hubble's observation and its interpretation in the framework of General 
Relativity (GR). This was a spectacular break through in our understanding (or misun- 
derstanding) of the Universe, and it became a major task to Cosmology and Astrophysics 
to explain this unexpected fact. Staying in the domain of classical GR and, consequently, 
considering the Friedmann's equations as valid, the only way to explain such present accel- 



the Universe is by considering the existence of some negative pressure dark energy 



eration o 

aai 

p|.This cosmic dark energy opposes the self-atraction of matter and is causing the 
expansion of the universe to be positively accelerated [4]. The most obvious candidate to be 
such dark energy is the cosmological constant and/or the vacuum quantum fluctuations of 
fields, which do have negative pressure. However, theorists estimates that the zero point en- 
ergies of the quantum fields must be at least 55 orders of magnitude larger than the critical 
density value. Hence, there must exist some yet unknown profound theoretical reason for 
the many contributions to the effective value of the cosmological constant be cancelled out 
to yield a number 55 orders of magnitude less then expected, or even zero. This is known as 
the cosmological constant problem. Some theorists believe that some profound symmetry 
requirements can be found to explain an exact cancellation, but not a partial one with an 
extreme fine tuning. In the case where the effective cosmological constant is indeed exactly 
zero, there were proposed some candidates in order to explain the present accelerated ex- 
pansion of the Universe as, for example, a very light, evolving scalar field called quintescence 

Another way to tackle this problem is by considering that presently, at cosmological scales, 
classical GR is not valid. In other words, instead of changing the right-hand-side (RHS) of 
Einstein's equations by introducing some new negative pressure fluid, one could try to find 
physical reasons which justify the modification of its left-hand-side (LHS) accordingly. How 
this can be done? 

In early works Q^l? a quantum minisuperspace model containing a free massless scalar 
field minimally coupled to gravity in a FLRW geometry was studied. These models were 
interpreted in the framework of the ontological Bohm-de Broglie (BdB) interpretation of 
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quantum mechanics, ljj] 11| Q], in order to extract predictions from the wave function of the 
Universe. This interpretation avoids many conceptual difficulties inherent to the application 
of the Copenhagen interpretation to the quantization of the whole Universe, where no place 
for a classical domain exists. The BdB interpretation does not need a classical domain 
outside the quantized system to generate the physical facts out of potentialities (the facts 
are there ab initio), and hence it can be applied to the Universe as a whole 1 . The solutions 
of the Wheeler-DeWitt equation for such scalar tensor model contain positive and negative 
frequency modes, the first leading to an expanding universe, and the second to a contracting 
one. There were constructed some particular superpositions mixing negative and positive 
frequency modes. In Ref. gaussian superpositions were studied and, for the case of flat 
spatial section, the Bohm guidance equations were reduced to a dynamical system. In this 
way the quantum trajectories were studied, emerging the following three kind of scenarios: 
periodic solutions, representing oscilating universes, bouncing universes, and models with a 
big bang followed by a big crunch. The bouncing universes contract classicaly from infinity 
until a minimum size, where quantum effects become important acting as a repulsive force 
avoiding the singularity, expanding afterwards to an infinite size, approaching the classical 
expansion as long as the scale factor increases. For the periodic solutions, the quantum 
effects are always important, and they do not grow enough to yield a large Universe as ours. 
The models with a big bang followed by a big crunch behave as the classical solutions for 
small values of the scale factor, but display quantum behaviour for large scale factor. These 
quantum effects are responsible for the turning over of these solutions from decelerated 
expansion to contraction. Near the big crunch, the quantum effects are again negligible. 
Bohmian trajectories which behave classically for small scale factors but quantically for 
large scale factors where already found in Ref. [jj. This is not surprising as it is well 
known that a large universe behaves classically or quantically depending on its initial 
quantum state. After these remarks, the natural question one can ask is if it is possible that 
quantum cosmological effects at large scales can mimic a negative pressure fluid and yield a 
positive acceleration for the whole Universe. The aim of this paper is to show with a simple 
model that it is indeed possible for some suitable initial quantum states of the universe. We 



1 Other alternative interpretations can be used in quantum cosmology, as the many worlds interpretation 
of quantum mechanics 
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take the flat model considered in Ref. j]J, and we consider another gaussian superposition 
of negative and positive modes solutions of the Wheeler-DeWitt equation. We write the 
Bohm guidance equations, which are reduced to a dynamical system, and we analyze the 
bohmian trajectories in configuration space. We find the two following scenarios, depending 
on the initial conditions: oscillating universes without singularities and with relative small 
amplitudes of oscillation, and universes which arise classically from a singularity, experience 
quantum effects in the middle of its expansion, and recover its classical behaviour for large 
values of the scale factor. We concentrate our attention on these solutions and we study the 
epoch where the quantum effects are important. We calculate its acceleration and explore 
its behaviour as a function of the scalar field <fi and of the logarithm of the scale factor, 
a = ln(a). We find that a positive acceleration of the universe can be obtained in such 
models, whose value can be adjusted by the choice of the free parameters of the model. 
This positive acceleration is a quantum effect. The mechanism is driven by the quantum 
potential, which appears in the modified quantum Einstein-Hamilton- Jacobi equation and 
modifies the usual classical trajectories. In this model, the acceleration is not forever: in the 
future, the universe recovers its classical deccelerated expansion. In this way, we present a 
possible alternative explanation for the accelerated expansion of the Universe today. Note 
that this explanation is based on quantum effects not only present in the scalar field, as 
described within a different approach in Ref.Q], but also in the geometry itself. 

The article is organized as follows. In Sec. II, we describe the classical model and we 
quantize it. In Sec. Ill, we introduce the Bohm-de Broglie interpretation of the quantized 
minisuperspace model presented in Sec. II. We study gaussian superpositions of the quantum 
solutions previously found, and we obtain the corresponding quantum bohmian trajectories. 
In Sec. IV we analyze the beahviour of the acceleration of the scale factor in the quantum 
bohmian trajectories, first qualitatively, by showing some period in the history of the model 
where the acceleration of its expansion is positive, then quantitatively, by comparing the 
curve relating the luminosity distance with redshift in the quantum model with the corre- 
sponding curve coming from the classical model suplemented by a cosmological constant. 
Sec. IV is for discussions and conclusions. 
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II. CLASSICAL AND QUANTUM MINISUPERSPACE MODELS 

In this section we make an overview of the models studied in Ref. [jj 

A. Classical Model. 

We start from the Lagrangian 

1 



L = y r ~g 

We consider the FLRW metric given (in isotropical coordinates) by 



(1) 



ds 2 



-N 2 dt 2 + 



a (ty< 



1 + | r 2)2 



dr 2 + r 2 (d0 2 + sin 2 



(2) 



the quantity e being the spatial curvature with values 0, 1, —1 for flat, spherical and hy- 
perbolic spatial sections, respectively. This line element will give, after inserting it in the 
Lagrangian (JJJ, the following action: 



S 



3V f Na 3 



4tt1 2 p 



-a 2 1 
+ 



2 \N 2 a 2 2GN 2 a< 
where we have set h = c = 1 and ' = 4. The quantity V is the volume divided by a 3 of the 
spacelike hypersurfaces, which are supposed to be closed, and lp is the Planck length. The 
total volume V depends on the value of e and on the topology of the hypersurfaces. For e = 0, 
V can have any value because the fundamental polyhedra of closed e = hypersurfaces can 
have arbitrary size For the case e = 1 and topology S 3 we have V = 2tt 2 . Defining 

2 Airll 



+ 



dt. 



(3) 



4 

3V ' 



'12, and omiting the bars, we obtain for the Hamiltonian: 



H = N 



P 2 f + P 2 ^ 
2a 



i 

2a 3 



(4) 



where p a and ps are the moments canonically conjugate to a and respectively, given by: 



Pa 



aa 
a 3 4> 



(5) 



(6) 
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A dimensionless scale factor is defined by a = ^ and the Hamiltonian becomes, omitting 
the tilde, 



2/3 



ea 



a 



a;' 



(7) 



As (3 is a multiplicative constant in the hamiltonian, we can set f3 — 1 without any loss of 
generality, keeping in mind that the scale factor which appears in the metric is a p h ys = (3a, 
not a. Defining now a = ln(a), we simplify the Hamiltonian obtaining: 



H 



N 

2^ 



(8) 



where 



Pa = 



e 3a a 



(9) 



P<t> 



e Aa <p 



(10) 



This Hamiltonian does not depend explicitly on <p. Hence, is a constant of motion, which 
we will call k. The classical solutions in the gauge N = 1 can now be listed: 



1. Flat model, e = 0. 

In configuration space, the classical solutions are: 

(j) = ±a + c 1 , (11) 
where C\ is an integration constant. In terms of cosmic time t they read: 

a = e a = (3fct)3, (12) 

0=^ + c, (13) 

These are solutions forever contracting or expanding from a singularity, depending on the 
signal of k, without any inflationary epoch. 
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2. Spherical model, e = 1. 



In this case we have, 



\ (14) 

cosh(20 - Cl ) ' K J 



where c\ is an integration constant. Conservation of p$ implies that 



k = e 3 >. (15) 

These solutions describe universes expanding from a singularity till a maximum size and 
contracting again to a big crunch. Near the singularity, these solutions behave as in the flat 
case. There is no inflation. 



3. Hyperbolic model, e = — 1. 



The classical solutions in configuration space are: 



a = e a = , — - -r, (16) 

|sinh(20-d)|' v ; 

where C\ is an integration constant. Again, from the conservation of ps we get 



k = e 3a <p. (17) 

These solutions describe universes contracting forever to, or expanding forever from, a sin- 
gularity. Near the singularity, these solutions behave as in the flat case. There is no inflation 
phase. The cosmic time dependence is complicated in the cases 2 and 3 and we will not 
write it here. 



B. Quantization. 

Let us quantize the model following the Dirac procedure ^|. The constraints become 
conditions imposed on the possible states of the quantum system. The operator version of 
the Hamiltonian (jSJ), obtained by setting a — > — i-jr^ and 4> — > — %-§r, must annihilate the wave 



S 



function Choosing a factor ordering which make it covariant through field redefinitions, 
the quantum constraint, i.e. the Wheeler-DeWitt equation, reads 



a^ + W + ee ^ ' (18) 



whose general solution can be written as 

/oo 
F{k)A k {a)B k {4>)dk, (19) 
-oo 

where k is a separation constant which in the classical limit corresponds to k, F(k) is an 
arbitrary function of k, the function B k (<p) reads 

B k (<p)=b 1 e ik * + b 2 e~ ik \ (20) 
and, for e = 0, the function A k (a) is given by 

A k (a) = a ie ika + a 2 e~ ika , (21) 

while for e = 1 it is 

e 2a e 2a 
A k [a) = aJik/2— + a 2 K ik/2 —, (22) 

and for e = — 1 it reads 

e 2a e 2a 
A k (a) = axJik/2— + a 2 N ik / 2 —. (23) 

The functions J, N, I, K are Bessel and modified Bessel functions of first and second kind. 

III. THE BOHM-DE BROGLIE INTERPRETATION OF THE QUANTUM 
MODEL 



The Bohm-de Broglie 
be summarized as follows 



BdB) interpretation of homogeneous minisuperspace models can 
1^ | : the Wheeler-DeWitt equation is 

H\p a {t),q a {t))*(Q) = 0- (24) 
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Substituting the wave function in polar form, \1/ = Ae tS , we have a complex equation and 
its real part produce, after dividing it by A: 

2 UM H + U{%) + Q{%) = °' (25) 
where Q{q,j) is the quantum potential, given by 

«*> = -2l^aSb' (26) 

In the BdB approach, the trajectories q a (t) are supposed to be real, independent of any 
observations. Equation is the Hamilton- Jacobi equation for them, but with an extra 
term given by the quantum potential Then we define 

P a = (27) 
dq Q 

where the momenta are related to the velocities in the usual way: 

,r = (28) 

In order to obtain the quantum trajectories we have to solve the Bohm guidance relations, 
P = IT ' wn i cn are ; i n this case, given by: 

Equations (|29|) are invariant under time reparametrization. Hence, even at the quantum 
level, different choices of N(t) yield the same spacetime geometry for a given nonclassical 
solution q a (t). There is no problem of time in BdB interpretation of minisuperspace quantum 
cosmology 2 . 

In the case of hamiltonian (|8j) . we obtain for the guidance relations (J2*9*|) 

dS *"* (30) 



da N 

OS e 3a d 



d(j) N 



(31) 



2 This is not true for the full superspace, see 0], although the theory remain consistent, see [i^, [2l| . 
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The modified Hamilton- Jacob i equation ()25|) reduces to 



> J 



da) 



-2 e 



la 



1 /d 2 A d 2 A 
~ 2A\'dd^ ~ da 2 



(32) 



where the last term in the LHS represents the quantum potential 



Q{a, 



1 

2A 



d 2 A d 2 A 



da 2 



(33) 



We will now apply the BdB interpretation to our minisuperspace model. We will restrict 
ourselves to the case of flat spatial sections, hypersurfaces with e = 0. In Ref. if, the 
following Gaussian superpositions of the solutions (|T9|) were studied: 



*x(a,0) = / F{k)B k {<t>)[A k {ot) + A_ k (a)}dk, 



(34) 



and 



* 2 (a,0) = J F(k)A k (a){B k (<f)) + B_ k (<p)]dk, (35) 
both with a2 = 62 = in Eqs ()20l21j) . and where the arbitrary function F(k) is the Gaussian 



F(k) = exp 



(k-d) 



2i 



a* 



(36) 



While in paper Q the solution ^ was studied in detail, here we will concentrate our analysis 
on the solution ^! 2 - Integrating (J5Sj) in k we obtain, for 



(a+(f>) 2 v 2 
4 



\&2 = aibi I o I v^rj exp 

exjp(id(a — </>)) 



exj>(id(a + </>)) 



+ exp 



(a-<j>) 2 a 2 
4 



(37) 



To obtain the quantum trajectories, we have to calculate the phase S of the above wave 
function and substitute it into the guidance equations (J30|) and (j3Tf . We will work in the 
gauge N — 1. Computing the phase of ^2, we obtain S = da + arctan( 2L y^) tan(— d<fi) 
which, after substitution in Eqs (j3()|31j) . yields a planar system given by: 



a 



So 2 sin (2d(f)) — 2dcos(2d(f)) — 2c?cosh(cr 2 a;^ 
e 3a 2[cos(2rf0) + cosh(a 2 a0)] 



(38) 
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• aa 2 sm(2d<j)) + 2dsmh(a 2 a<j)) 

^ = ~~ e 3 «2[cos(2ci0) + cosh(a 2 a0)] ^ 
Equations ()38I39|) give the direction of the geometrical tangents to the trajectories which 
solves this planar system. By plotting the tangent direction field, it is possible to obtain the 
trajectories. The vertical line = divides the configuration space in two symmetric regions. 
The line a = contains all singular points of this system, which are nodes and centers. The 
nodes appear when the denominator of the above equations, which is proportional to the 
norm of the wave equation, is zero. No trajectory can pass through these points. They 
happen when a = and (f> = (2n + 1)^, n an integer, with periodicity 7c/\d\. The center 
points appear when the numerators are zero. They are given by a = and = |f cot(<i0). 
As |0| — > oo these points tend to rnr/d (zeros of tan(<i0)). As one can see from the above 
system, the classical solutions (a(t) oc t 1//3 ) are recovered when — > oo or \a\ — > oo, the 
other being different from zero. 

We present a field plot of this planar system in Fig. Q for the case d — — 1, a — 1. 
Depending on the initial conditions, we can see two different possibilities. Near the center 
points there are oscillating universes without singularities and with amplitudes of oscillation 
of order 1. The other possibility is given by non-oscillating universes. A non oscillating 
universe arises classically from a singularity, experiences quantum effects in the middle of 
its expansion, and recover its classical behaviour for large values of a. 

IV. THE ACCELERATED EXPANSION. 
A. Qualitative approach 

In this subsection we show how our approach for the explanation of the accelerated 
expansion works in a qualitative manner, i.e., with the parameters of the wave packet (d 
and a) adapted for readable numerical treatment, without any fitting with usual cosmo logical 
orders of magnitude. The aim is just to show that in some period in the history of such 
quantum models where the expansion is positively accelerated. We take here d = — 1, 
a = l for the numerical computations. The quantum effects appearing in the middle of 
the non periodic bohmian trajectories described above can deviate them from their classical 
decelerated expansion to an accelerated one. We will show that this is indeed the case of 
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FIG. 1: Field plot of the system of planar equations (|38|39jl coming from the wave function (|37|) .For 
numerical simplicity we choose the values —d = a = 1. We see the two possibilities: trajectories 
corresponding to oscillating universes without singularities and trajectories corresponding to non- 
oscillating universes that come from a singularity, experience quantum effects in the middle of their 
expansions, and recover their classical behaviours for large values of a. 

this model. From a = lna(t) we have 



From Eq. (jSHj), oc can be viewed as a function of the canonical variables at = f(a, <p). Then 
we have 




(40) 



a 




Computing the derivatives with respect to a and 0, and substituting at and (j> from equations 
(|HH|) and (}3T?|) . respectively, we obtain 
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(cx 6 3 + 4(x 2 d 2 - a 6 (j)a 2 ) (sin(2d0)) 2 sinh(a 2 a0) + (2a 4 d0a - 8(x 2 d0) sin(2d0) (cos(2d0)) 2 + 
(2a 4 d0a - 16(x 2 d0) sin(2d0) cos(2d0) cosh(a 2 a0) + 24d 2 (cos(2d0)) 2 cosh(a 2 a0) + 
24d 2 cos(2d0) (cosh(a 2 a0)) 2 - 80d<x 2 sin(2d0) (cosh(a 2 a0)) 2 + 4<x 2 0d 2 (cos(2d0)) 2 sinh(a 2 a0) + 
4cr 2 0d 2 cos(2d0) cosh(o" 2 a0) sinh(cr 2 a0) + 2<r 4 (sin(2d0)) 3 da — 2<7 4 0sin(2d0) (sinh(<7 2 a:0)) ad + 

8d 2 (cos(2d0)) 3 + 8d 2 (cosh(a 2 a0)) 3 (e 3a ) 2 (cos(2d0) + cosh(a 2 a0) 



"(42) 



The equation above gives the acceleration - as a function of a and 0. If one integrates 
the system ()38I39|) to obtain = 0(a), the quantum version of the classical equation (i.e. 
Raychaudhuri equation for the Friedmann model) 

a AtiG , „ « . • . 9 1 , 

- = + 3 P « -W> « —6 43 

a 3 a b 



can be obtained. Taking the limit where the absolute value of a is very large in Eq. 
one recovers the classical behaviour given in Eq. (}4~3*|) . - oc — 1/a 6 . Note that the quantum 
analog of the Friedmann's equation can also be easily obtained from Eqs. ()30l31l32|) . yielding 
(recovering the unities) H 2 = (ft 2 + j3 4 c 2 Q(a, 0)/ap hys . 

We can represent - in a tridimensional plot as a function of a and 0. In this plot we can 
see the regions on the plane a — in which the acceleration is possitive, negative or zero. 
We show this plot, for the parameters d — — 1 and o = 1, in Fig. El One can see the classical 
behaviour a/a oc —1/a 6 for a — > (a — > — oo) and a — > oo (a — > oo), but near the region 
a = 1 (a = 0), a clear departure from classical behaviour is observed, and positive values of 
a/a are obtained. Fig. El shows the acceleration a/a as a function of a for = 1.66. 

As we pointed above, the quantum potential is the cause of this positive acceleration. 
Fig. El shows the quantum potential in the a — plane where we can see that it is different 
from zero in the regions of positive acceleration. A trajectory passing through this region on 
the plane a — will correspond to a universe experiencing an accelerated expansion. Fig. El 
shows the quantum potential as a function of a for = 1.66. Note also the increasing in the 
acceleration and in the quantum potential as we approach the node point a = 0, = tt/2. 
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1 



FIG. 2: Acceleration a/a as a function of 4> and a. For numerical simplicity we choose the values 
-d = a = l. 
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FIG. 3: Acceleration d/a as a function of a for <fi = 1.66. For numerical simplicity we choose the 
values — d = a = 1. 
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FIG. 4: The quantum potential as a function of <p and a. We see that it is different from zero in 
the region of positive acceleration. For numerical simplicity we choose the values — d = a = 1. 




FIG. 5: The quantum potential as a function of a for <j) = 1.66. We see that it is different from zero 
in the region of positive acceleration. For numerical simplicity we choose the values —d = a = 1. 
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B. Comparison between the quantum acceleration of the Universe's expansion 
and the presence of a cosmological constant 



We will now compare quantitatively the accelerated expansion of the Universe caused by 
such a quantum cosmological effect with the one generated by the presence of a cosmological 
constant in the classical model . First of all, we must recover the units in Eqs. (j38l39j) . 
This is done by multiplying the RHS of these equations by c//3, where j3 = y~^lp and 
V = Vph ys (t) / ap hys (t) , where a p h ys = o,f3 is the physical scale factor and Vphys is the physical 
volume of the spacelike section. Then we obtain 

H = a = f(a,<f>) VP , (44) 

y p hys{t)tp 



where 



g{a,4>) * , (45) 



t( m\ 0a 2 sin(2#) 

2|cos(2a</>) + cosh(a z a(p)\ 

aa 2 sm(2d<ft) + 2d sinh(<T 2 a0) . . 

9 ^ a ^)= 2[cos(2d0) + cosh( ( x 2 a0)] ' ( ' 

Vp = ^lp is the Planck volume, and tp is the Planck time. We will compare this quantum 

cosmological model with the original classical free scalar field model, classically equivalent 

to stiff matter, with flat spatial section, suplemented with a cosmological constant as an 

alternative source for accelerated expansion. This classical model satisfies the Friedmanns 

equation 

H 2 = + A = H%[{1 - 00(1 + zf + Ov], (48) 

where is a constant such that the energy density of the field is = c<p/a 6 , Q\ = 
A/(8vrGp crit ), p cr it is the critical density, and Ho is the Hubble's paremeter today. The 
deceleration parameter today, q = — a/(aHo), is given by 

g = 2 - 3fi A . (49) 
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To obtain the luminosity distance as a function of z one can calculate numerically the integral 



d L = (1 + z) / — - = — - — / (50) 
Jo H(y) H Jo J(l-0 A )(l +2/ 6 + o A 



As a power series, it can be written as 

H d L = z + 2 ?(-l + ^n A ^j + z 3 (l - fi A ) Q - + (51) 

In the quantum cosmological problem, we have to deal with Eqs. (|44|45|) . There are 
four arbitrary parameters: a, d and two coordinates in the (a, <$) plane, 0o designating 
the particular trajectory in Fig.(£Q) which represent our Universe, and «o designating the 
present moment in a particular trajectory. The present value of the Hubble's paremeter in 
a particular trajectory coming from Eq. f)44j) is given by 

Ho = f(a ,<p )-^- (52) 

where 

ft a\- o a 2 sin(2rf0 o ) 

J\OLo,(po) = — — - — -tt-« (53) 

2[cos(2a</> ) + cosh(cHa o 0o)J 

In order to obtain a model similar to our present Universe whith Hq — 10 _18 s _1 and Vp° hys > 
10 82 cm 3 , one must have 

ft m\- O ^ 2 sin(2c/0 o ) , in i2o fKA \ 

f(a ,(j) ) = — 10 54 

2(cos(2rf0 o ) + cosh(a 2 a;o0o)) 

This huge number can be obtained by choosing a very large value for a 2 (the gaussian 
in the wave function would be almost flat indicating no preference in the choice of k, or 
physically, no prefered choice in the strength of the initial explosion), a very large value for 
\d\ (a gaussian centered in a very negative value of k, or a choice for a very strong initial 
explosion), or trajectories passing very close to the node point, where the denominator of 
the above expression approaches zero. 

In the case of a large value of a 2 , choosing «o — 1/c 4 — 0, one can check from Eq. (j42|) 
that 

a. 1 2 cos(2<i0o) + 1 , . 

qo = ~a l °H 2 ^ cos(2# ) + l 1 } 
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. For 2d0 o = 2nn + x with x G (2.145, 2.15) one has -0.21 < q < -0.17. Note that, for the 
classical model with Q A = 0.73 22], q = -0.19 (see Eq. flU). 

The supernovae measurements relate the luminosity distance dz with z. Hence, it would 
be instructive to compare the quantum cosmological luminosity distance d\{z) 

Jo f(a(y),<t>(y)) c 

with the one given in Eqs. (|50|51jl . We will expande G(y) = e 3a( - y ^ / f(a(y),cj)(y)) in powers 
of y around y = (today) 



dG y 2 d 2 G y 3 d 3 G 

&y~ ly=0 + 2! d7 ly=0+ 3! d^ 

up to third order. The operator d/dy can be written as 



G(y) = G(0) + y— 1, =0 + ^ — U + &"37s ly=o + - (57) 



dy (l + y) da (l+j/)U da^/' 1 ; 

where d0/da can be easily obtained from Eqs. ()38J39j) and a = ln(a) = — ln(l + y) + const.. 
Recalling that 

HoPaf) _ H a 3 hys _ -^oQph ys ^p3V _ H tpV phys ^ 1 nl20 

C C/3 2 47T Vp 

we obtain 



10 120 , (59) 



H d q L 



10 



120 



/o 



: . 2 fl + A^lo 



2ag d;z 



3 / /o dG /o d 2 G 
+ z ( ^—5- -r- 1 + — ^ 77 1 



V2ao dz 6a 3 , dz 2 







+ 



,, /o d 2 G /o d 3 G , 
^ ,_ ^ |o + 24aId^ |o) + -" 



(60) 



,6a 3 , dz 2 

Comparing Eq. (j60|) with Eq.f|51jl. one can see that /o must be of order 10 120 , as we have 
already concluded in Eq. (j53jh Here we will choose \d\ as the big number to make fo ~ 10 120 , 



and |«o| — 7/ Ml? 7 being an arbitrary number of order 1. This is to assure that 4Jr|o is of 
order 1 or less for any n > 1, which is not the case if we choose cr ^> l 3 . Using Eq. (|58j) . we 
can calculate the coefficients in Eq.(|5U|l. obtaining 



3 In this case, the series expansion of Ea. 1)56(1 is not meaningful at lower orders, and the integral must be 
performed by other methods. 
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z 



FIG. 6: The luminosity distance as a function of redshift. The thin line curve corresponds to the 
quantum model, the dotted curve is for the classical model with a cosmological constant, and the 
thick line curve is for the classical model without cosmological constant. 



2 12 6[1 + cos(2d0 OJ 



, z 4f 1 0|]a 6 7sm(2#o) \ 

12 3[l + cos(2c/0 o )] 2 2[l + cos(2c/0 o )] 4 i 

In order to obtain Eq.([61jh we have used that \d\ » 1 and |o? | = j/\d\ is very small 4 . 

There are many values of <po, a and 7 which makes the graphic of H Q d q L similar to the 

one obtained from Eq. (jH0~j) for = 0.73. For instance, for 7 = 1 we set 2d(f) = 0.64 ± 2mr 



4 Note that the first correction to the classical model without cosmological constant (see Ea. (|51ll with 
Clji = 0) come in the cubic term. The correction to the deceleration parameter at this particular moment 
is negligible with this choice of parameters. However, the corrections in the cubic and forth terms can 
be adjusted in order to make the curve obtained from Ea. H61|) close to the corresponding curve obtained 
from Ea. H50l) . as we will see. 
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and the coefficients of z 3 and z A become, respectively, 2.6 and —2.25. 

In Fig® we show a plot of HocIl(z) given by Eq.(|5Uj). H cIl{z) given by Eq.(J3UJ) with 
Q Lambda = 0, and H d q L (z) given by Eq. fJBT]) with the coefficients of z 3 and z A being equal 
to 2.6 and —2.25, respectively. Note that for small values of z they are close but, for 
intermediary values of z, the quantum d\{z) remain close to the cosmological constant 
cIl(z) while both separates of the pure stiff matter di(z). Of course, for bigger values of 
z, the quantum d 9 L (z) may separate strongly from the cosmological constant di{z). Hence, 
quantum comological effects may mimic a cosmological constant in some region but not 
everywhere. The two models are distinguishable. 

V. CONCLUSION 

We have studied gaussian superpositions of positive and negative frequency mode solu- 
tions of the Wheeler-DeWitt equation corresponding to a scalar-tensor model in minisuper- 
space in the case of flat spatial section. According to the Bohm-de Broglie approach to 
quantum cosmology, the quantum trajectories representing dynamical universes evolving in 
time were studied. We have shown that it is possible to have universes which arise classi- 
cally from a singularity, undergo a positive acceleration in the middle of its expansion, and 
recover its classical behaviour for large values of the scale factor. We have shown that this 
positive acceleration, which can be made compatible with observations for many choices of 
initial conditions, is due to a quantum cosmological effect driven by the quantum potential, 
according to the BdB interpretation of quantum cosmology. In this way, it may be possible 
to explain the positive acceleration suggested by the recent measurements of high redshift 
supernovae without postulating a new contribution to the energy density of the Universe as 
the dark energy. Note that this acceleration is caused by quantum effects not only present 
in the scalar field, as described in Ref . , but also in the geometry itself. We consider the 
Universe as a quantum system no matter its size. It is possible to have small classical uni- 
verses and large quantum ones: it depends on the state functional and on initial conditions 

QQ. 

The quantum cosmological explanation for the acceleration of the Universe presented 
in this paper needs to be studied further, not only because it would be an alternative 
explanation for a misterious behaviour of the present Universe without appealing to any 
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new form of energy, but also because it is a possibility of an observable physical effect of 
quantum cosmology. Furthermore, quantum cosmological explanations may be suported by 
symmetry principles which are absent in the classical domain. As an example, we have 
seen that the huge cosmological numbers may be explained by some "democracy" principle 
stating that any value of the velocity of expansion (the constant k) is equally good (a is 
very big, the gaussian is almost flat). Of course, more elaborated models taking into account 
relevant matter sources like dust and radiation must be studied. This will be the subject of 
our future investigations. 
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